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ABSTRACT 

How did the dwarf galaxy population of present-day galaxy clusters form and grow over time? 
We address this question by analysing the history of dark matter subhaloes in the Millennium- 
II cosmological simulation. A semi-analytic model serves as the link to observations. We 
argue that a reasonable analogue to early morphological types or red-sequence dwarf galax- 
ies are those subhaloes that experienced strong mass loss, or alternatively those that have 
spent a long time in massive haloes. This approach reproduces well the observed morphology- 
distance relation of dwarf galaxies in the Virgo and Coma clusters, and thus provides insight 
into their history. Over their lifetime, present-day late types have experienced an amount of 
environmental influence similar to what the progenitors of dwarf ellipticals had already ex- 
perienced at redshifts above two. Therefore, dwarf ellipticals are more likely to be a result of 
early and continuous environmental influence in group- and cluster-size haloes, rather than a 
recent transformation product. The observed morphological sequences of late-type and early- 
type galaxies have developed in parallel, not consecutively. Consequently, the characteristics 
of today's late-type galaxies are not necessarily representative for the progenitors of today's 
dwarf ellipticals. Studies aiming to reproduce the present-day dwarf population thus need to 
start at early epochs, model the influence of various environments, and also take into account 
the evolution of the environments themselves. 

Key words: galaxies: abundances - galaxies: clusters: general - galaxies: dwarf - galaxies: 
evolution - galaxies: haloes - galaxies: statistics 



1 INTRODUCTION 



, The filamentary large-scale distribution of galaxies taott et al] 
■ l2005i) implies that any given galaxy cluster is connected to a 
number of surrounding environments with different characteristics. 
, These connections evolve with time: galaxies and galaxy groups 
get accreted, c ontributing to the growth of cluster and supercluster 
environments ( Springel et al ; Klvpi n et al 1). As a conse- 
quence of this "cosmic web" and its inhomogeneous density distri- 
bution on megaparsec scales, observed galaxy properties at low red- 
shift depend not only on the local env irormient - their current group 
or cluster (e.g. lGavazzietaT1l2010l) - but also on the surround- 
ing la rger environment ( Mahaian et al.ll2oTll 1201 2l : iLietzen et al.l 

To investigate the influence of environment on galax- 
ies, dwarf galaxie s serve as id eal tracers. Their abundance 
is large (T rentham & Tullvl l2002h . and their shallow potential 
jde Blok et aU2 008) makes them susceptible to ex ternal processes , 
namely gravitational forces and ram pressure (iGnedinI l2003bl : 
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|Mayer et al."200^. However, observations can merely provide a 
present-day snapshot of different dwarf galaxy populations that 
have evolved in parallel, not a sequence of their past evolution. We 
therefore utilise cosmological simulations and models to gain in- 
sight into the history of dwarf galaxy populations within a ACDM 
framework. Our focus lies on present-day massive galaxy clusters, 
which are the largest and probably best studied enviro nments of the 
nearb y u niverse beyond the L ocal Group (see, e.g., iBoselli et al.l 
I2OIII and lFerrarese et al.ll2012h . 

Observations of dwarf galaxie s in nearby clusters have 
become increasingly refi ned dAdami e t al. 2005; Ferrarese et al, 
2006; McDonald et alj boogi; IPennv et alj 2009; Hamm er etal. 
201a ; .Liederet al.„2012l ; lRvs et al.ll2012l ; ISmith et al...2012^ . and 
by dedicated sirnulations and m od 



hav e been comp lemented 

els (Moore et al. 'l999'; 'Frugeen & De Lucia" '2008'; Boselli et 
i2008b; Aguerri & Gonzalez-Garci'a 2009; De Rijcke eta l. 201 
ISmith et al.ll20ld ; ISchroven et al .1120 111 : I Smith et al.H2012bl) . How- 
ever, ACDM N-body simulations of cosmological volumes have 
only recent ly reached the regime of dwarf galaxies in their mass 
resolution teovlan-Kolchin et alj |2009| ; iKlvpin et alj 1201 ih . The 
large dynamic range between a dwarf galaxy and its host cluster 
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makes cosmological hydrodynamical simulations of their evolution 
difficult. A viable alternative is provided by semi-analytic models 
applied to merger trees of N-body simulations. While such mod- 
els have contributed to increase our understanding of higher mass 
galaxies, they have not been exp loited much in the study of clus- 
ter dw arf galaxies until recently dPeng et al.ll2008t iDe Lucia et al.l 
l2012 b). paralleling a larger number of studies for Milky Way sized 
haloe s (e.g. Maccio et al.ll20ld iFont et al]|201lL and Helmi et al.i 
l2012h . The first semi-analytic model b ased on the Mil lermium-II 
simulation ( iBovlan-Kolchin et alj200sl) was released bvl Guo et al.l 
( I2OI ih . It was tuned to reproduce the redshift-zero stellar mass 
function down to log(M*/M0) ~ 7.5 and the luminosity func- 
tion dow n to Mr = —15 mag. In this series of papers, we use the 
model of lGuoetalJ ( l201lh in combination with observational data 
for nearby galaxy clusters to investigate the properties of their faint, 
low-mass galaxy population. 

Of particular interest are early-type dwarf galaxies, as they 
are the d ominant galaxy populatio n by number in massive galaxy 
clusters djerien & Tarmnanrilll997h . Their origin is still not fully 
understood, for several reasons. First, it is still debated whether or 
not their basic scaling relations are a mere continuation of those 
of bri ght earl y types or ar e distinct from the m (Jerjen & Bin ggeh 



1997 : Graham & Gu zman 20 031: ICote et al.l 120071: IJanz & Liskei 
200g ; IChilingarianI i2009.; Janz & Lisked l2009l: iKormendv et al 



2009l : lGlass et al.l20Ill : lGrahaml201ll : lKormendv & Bendej2012h . 
Second, despite their similar overall characteristics in different 
environments (De Rijcke et al. 2009), an increasing number of 
studies has shown a pronounced complexity in the structure, 
internal dynamics, and stellar populations of early-type dwarf 
galaxies, revising th e pictu re of a sir nple galaxy population 
(Binggeli & Cameron 1991; Geha et al. 2003; Michielsen et al. 
L2OO8; Tully_&_Tren tham 2008; Lisker etal. 2 Q09; Toloba et al. 
I2OI1I ; iKoleva et alj i20Ill ; IPaudel et alj I2OI1I ; Ijanz et alj I2OI2I) . 
Third, due to their low luminosity, mostly faint surface bright- 
ness, and small size, studies of dwarf galaxy populations are 
mainly restric ted to the nearby universe (or at least to low redshift, 
lBarazzaetai]|2009i') . thus not probing a large range of lookback 
time. Fourth, it needs to be appreciated that the regime of poten- 
tial late-type progenitors of early-type dwarfs is also very diverse, 
rang ing from thin late-ty pe spiral galaxies with small or no bulge 
(e.g. iKautsch et ai]l2006l') , over blue compac t dwarf galaxies with 
starburst activity (e.g. IPapaderos et al.ll 1996b. to diffuse irregulars 



with low-level star formation (e.g. 



I Zeel200lh . All of these types 
coexist in the luminosity range of about 10* to 5 • 10® L© (e.g. 
ISandage & Bin ggeli 1984), which roughly corresponds to an abso- 
lute magnitude range of —15 < Mr < —19 mag. Fifth, late-type 
galaxies that we see today are those that have not been transformed 
to early types, and may thus be only partly rep resentative for the ac - 
tual progenitor population of early types (e.g. iBoselli et alj2008ah . 

The last point provides one of the main motivations for the 
current study. A number of environmental transformation processes 
have been investigated, whose combined effect may lead from ro- 
tating, gas-rich, star-form ing low-mass galaxies to dynamically hot , 
gas-poor, quiescent o nes dMoore et al.ll996t|Mori & Burkertl20"ool ; 
Ivan Zee et al] 12004 j; Ivon der Linden et al.ll2010h . Yet at which 
epochs do these have to operate, in order to reproduce the real 
present-day early-type dwarf population? What were the 'input 
galaxies' like at those epochs, and what were the characteristics 
of the environment they entered? To what extent can a galaxy's lo- 
cation at present tell us about the strength and duration of the envi- 
ronmental influence it experienced in the past? These questions set 



a basic framework - within a ACDM universe ~ for understanding 
the origin of the dwarf galaxy population that we observe today. 

In Weinmann et al. ( 201 1 , hereafter Paper I) we found that the 
dwarf galaxy abundances, velocity dispersions, and number density 
profiles observed in nearby massive clusters are generally well re- 
produced in the Millennium-II simulation and the lGuo et al. I boiih 
semi-analytic model. However, the comparatively low number of 
faint galaxies in the Virgo cluster center, within 300 kpc around 
the central galaxy M87, is not reproduced in any model cluster. 
The model may underestimate tidal disruption for faint galaxies, 
since the dwarf-to-giant ratio - defined in terms of luminosity - 
is systematically higher in the model than in the observed nearby 
clusters. We also found indications that the model probably overes- 
timates environmental effects that lead to star formation quenching 
in galaxy groups. 

In the study presented here, we focus on the mass loss and in- 
fall history of subhaloes in the massive clusters of the Millennium- 
II simulation. Semi-analytic model quantities are used mainly for 
selecting subhalo samples by galaxy magnitude, and for track- 
ing galaxies with tidally stripped subhaloes. As in Pape r I and 
in Guo et al. (2011), we use a WMAPl cosmology ( Sperg eTet al.l 
I2OO3I) and assume h = 0.73 throughout. 

The paper is organised as follows. Section |2] characterises 
the observational galaxy samples. The dark matter simulation and 
semi-analytic model are described in Section|3] Our analysis is pre- 
sented in Section |4] including the comparison between simulated 
subhalo populations and observed galaxy populations. This is fol- 
lowed by a discussion in Section|5]and by our conclusions in Sec- 
tion|6] 



2 OBSERVATIONAL SAMPLES 

Details of our Virgo, Coma, and Perseus cluster samples are pro- 
vided in Paper I. Here we briefly describe the sample selection and 
characteristics. All samples are limited to "dwarf" magnitudes by 
requiring the r-band absolute magnitude to be Mr > —19.0 mag. 
At the faint end, magnitude completeness limits were chosen to 
avoid losing galaxies with very low surface brightness. As outlined 
in Paper I, very compact galaxies may be missed in the Coma and 
Perseus cluster samples, but are not expected to contribute more 
than a few percent to the populat ion. All photometric values were 
corrected for Galactic extinction (ISchlegeletal.ll998b . 



2.1 Virgo Cluster 

Our Virgo sample is based on all certain a nd possible members 
listed in the Virgo Cluster Catalog (VCC, Bingge li et"!!] Il985l 
I1993I) . with membership updated by one of us (TL) in May 2008 
through new velocities given by the NASA/IPAC Extragalactic 
Database (NED), many of which were pro vided by the S loan Dig- 
ital Sky Survey (SDSS, Adelman-McCart hv et al.ll2007h . We ex- 
clude galaxies that are likely members of t he so-called M and 
W clouds, located at a distance of 32Mpc 1 Gavazzi et alj|l99^ 
in the western part of Virgo. Following iGavazzi et al. , galaxies 
in the projected region of these clouds and with a velocity wlg 
relative to the Local Group larger than 1900 km/s (with «lg ~ 
f'hciioccn. + 220 km/s) are assumed to belong to the clouds, and are 
therefore excluded from our sample. Note that this had not been 
done in Paper I, but only affects 34 of over 500 objects. Galax- 
ies for which no velocities are available remain in o ur sample. 
We use a distance modulus of m — Al — 31.09 mag jMei et al] 
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l2007l : lBlakeslee et al.ll2009l) for all sample galaxies, corresponding 
to d = 16.5 Mpc. With the adopted WMAPl cosmology, this leads 
to a spatial scale of 79 pc/" or 0.286 Mpc/°. 

Total r-band magnitudes and col o urs from uqriz-bwd& 
wer e measured by iLisker et all (l2007h . Ijanz & LiskeJ (|2009|) , 
and Meyer etal.l (l2013h on SP SS images from data release 5 
jAdelman-McCarthy et al.ll2007|). This inc luded a proper sky sub- 
traction method (described in Lisker et a l? 2007) that avoids the se- 
rious oyerestimation of the sky by the SDSS pipeline for nearby 
galaxies of large apparent size. For a small fraction of the sam- 
ple, r-magnitudes were obtained by transforming the VCC B- 
magnitudes (see the appendix of Paper I). The r-band completeness 
limit was estimated to be Mr < —15.2 mag, which we adopt as the 
limit for our sample selection. 

Our final working sample contains 521 galaxies, of which 442 
have spectroscopic heliocentric velocities. The sample is nearly 
complete out to a projected clustercentric distance of 1.5 Mpc from 
the central galaxy M 87 - comprising 404 galaxies - but includes 
galaxies up to 3 MpcQ 

To compute the red galaxy fraction, we split galaxies using the 
g — r colour cut from Paper I: 



{g - r)cut = 0.4 - 0.03 ■ (13 + Mr) 



(1) 



This cut was chosen blueward of the bulk of galaxies that populate 
the red sequence. 



2.2 Coma Cluster 

O ur Coma cluster sample is based on the SDSS data release 
7 jAbazaiian et al. I|2009). Total r-band magnitudes and ugriz- 
colours are extracted from the SDSS as Petrosian magnitudes. For 
all galaxies we us e a distance modulus of m — Af = 35.00 mag 
dCarter et aklbOOSh . corresponding to d = 100.0 Mpc. We initially 
include all galaxies with —19.0 < Mr < — 16.7mag (see Paper I) 
and out to a projected distance of 4.2 Mpc from the cluster cen- 
ter, which is defined as midway between the two massive central 
ellipticals NGC 4874 and NGC 4889. 

Objects for which SDSS redshifts are available with a redshift 
confidence of 95% or higher are considered as members if they 
lie in the range 4000 < cz ^ lOOOOkm/s. The SDSS spectro- 
scopic coverage reaches only to Mr < —17.3 mag. Those poten- 
tial member galaxies for which the SDSS does not provide spectro- 
scopic redshifts of acceptable quality were visually inspected to ex- 
clude contaminants (stars, parts of galaxies, and objects artificially 
brightened by a star's halo or a bright galaxy). The total number 
of objects is 1662, with 481 of them having reliable spectroscopic 
redshifts. 

We then correct the sample statistically for the number of con- 
taminating background galaxies. This is based on a radial number 
density profile (Appendix A3 of Paper I), which we inspect visually 
to adopt the value where the profile flattens out as our background 
value. We verify that this approach works rea sonably well by com- 
paring our selection to the galaxy catalogue of iMichard & AndreonI 
yOOe) for the Coma cluster center (see Paper I for a discussion). 

To estimate the red galaxy fraction, which will be used in Sec - 
tion 14.31 we first apply a k-correction dChilingarian et al.lbOld) 



^ In Paper I, we mistakenly stated that a restriction to < 1.5 Mpc 
would exclude the M 49 subcluster, while most of it is actually included 
)Binggelietai]ll987l) . However, none of the results from Paper I would 
have changed significantly if this subcluster would have been excluded. 



to the g ~ r colours of all galaxies, assuming their redshift is 
z = 0.023. The median k-correction is 0.036 mag. Then we 
background-correct red galaxies as described above, using the 
colour cut of eq.[T] 



2.3 Perseus Cluster 

Our Perseus sample is constructed in a similar way as the Coma 
sample, based on the SDSS data release 7 and using a statistical 
correction for background galaxies. SDSS spectroscopic redshifts 
are not available in this region. We use a distance modulus of rri — 
Ad = 34.29 mag for all galaxies, coiTesponding to a "Hubble flow 
distance" of d = 72.3 Mpc that is given by NED based on the 
heliocentric velocity of 5366 km/s from Stmble & Rood ( 1999). 

We initially include all galaxies with —19.0 < Mr < 
— 16.7 mag and out to a projected distance of 3.8 Mpc from the cen- 
tral galaxy NGC 1275. However, the SDSS coverage of the Perseus 
cluster outskirts is not complete. The incompleteness remains mod- 
erate out to 2.5 Mpc and is coiTected for each galaxy by a factor that 
takes into account the missing area at its clustercentric distance. 
This area is calculated for each galaxy's position as the part of an 
annulus around the cluster center that is not included in the roughly 
rectangular coverage of SDSS data. 

As for the Coma sample, all objects were visually inspected 
to exclude stars, artifacts, and parts of galaxies. The total number 
of galaxies in the sample is 1607. These are then corrected statisti- 
cally for the background level revealed by a radial number density 
profile (Appendix A4 of Paper I). While this level cannot be ac- 
curately determined and may influence the number counts in the 
cluster outskirts, its uncertainty is of minor relevance for the dense 
inner regions of the cluster. 

As for the Coma galaxies, we first apply a k-correction to 
the g — r colour before background-correcting red galaxies. With 
the adopte d redshift of z = 0.01 7, the median k-correction is 
0.029 mag dChiUngarian et al.l2010l) . 



3 DARK MATTER SIMULATION AND GALAXY MODEL 
3.1 Overview 

Semi-analytic models describe galaxy formation and evolution on 
the basis of simple analyti c recipes that are app lied to dark matter 
merger trees ( Kau ffmann et al .119931: ICole et al.ir200 0: Bower et al 



b006). In our study we use the semi-analytic model of IGuo et al 



( 1201 1,1 . which was applied t o the merger trees of the Mi llennium- 
II cosmological simulation jBovlan-Kolchin et alj|200sl ). The lat- 
ter simulates the evolution of the dark matter distribution within 
a periodic box of 137 Mpc side length and a particle mass of 
9.45 X 1 0*^ Mq. 

The Guo et al ('2011') model is based on the model of 
iDe Lu cia & Blaizot (2007), but contains several modifications. The 
supernova feedback efficiency for low-mass galaxies was increased 
considerably, in order to fit the stellar mass function of galaxies at 
the low-mass end. Furthermore, the prescriptions for calculating 
the sizes of galaxies and environmental effects were modified (see 
Paper I for details). In addition, the AGN feedback efficiency was 
increased. 

Friend-of-Friend (FoF) haloes were defined in the Millenium- 
II simulation by linking particles with separation below 0.2 of 
the mean value (Davis et al. 1985). Within these FoF haloes, the 
SUBFIND algorithm CSpringel et al..,2001i) identified subhaloes, to 
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which the semi-analytic model associated galaxies. Once the sub- 
halo mass falls below the mass of the galaxy, the model employs an 
analytic prescription for the galaxy's orbit (see Guo et al. 2011). 

The galaxies and galaxy clusters whose d ark matter haloes 
were simulated by iBovlan-Kolchin et alj j2009h and whose baiy- 
onic configuration was modeled by Iguo et al.l ( 1201 ih will be ref- 
ered to as "model galaxies" and "model clusters" throughout 
the paper. Their prope rties are publicly provided by the Virg o- 
Millennium Database iLemson & The Virgo Consortium! l2006h FI 
While we primarily make use of properties provided by the dark 
matter simulation itself (such as position and mass of a subhalo), 
we rely on the semi-analytic model for selecting galaxies by their 
r-band magnitude, their stellar mass, and for tracking so-called or- 
phan galaxies, which have already lost their dark matter subhalo. 
Magnitudes of model galaxies are based on the stellar evolutionary 
sy nthesis model s of Bruzual & Chariot (2003), which were found 
by iHansson et alj ( 120121 ) to well reproduce the ugriz-colors of 
galaxies in the nearby universe. Since there are indications that the 
semi-an alytic model does not c orrectly reproduce dust attenuation 
(Paper I: IWeinmann et alj|2012l) . we chose to use dust-free magni- 
tudes. 



3.2 Selection of model clusters and their galaxies 

For comparing the model galaxies and clusters to observations, 
we define three samples of model clusters. The full sample com- 
prises the 15 most massive galaxy clusters at redshift zerc[f| of 
the Millennium-II simulation. The "SAM-V" sample only uses 
the twelve least massive clusters of the full sample, thereby cov- 
ering the published range of virial masses of the Virgo cluster 
(1.4-4.0 X 10^"^ Mr based onlBohringer et al.ll994 iMcLaughlinI 



1 19991 ISchindleretai]|l999l lUrban et al.ll201 ll and Paper I). The 
"SAM-CP" sample only uses the three most massive clusters of the 
full sample, with 4.5 x 10^*, 4.7 x lO", and 9.3 x lO^'^ Mp. These 
cover the range of published Perseus cluster masses (lEvles et al.l 
I1991I : ISimionescu et aljboul ; Paper I), a nd reach at least close t o 
the Coma cluster mass of 1.3 X 10^^ Mq dLokas & Mamo3l2003h . 
When using observer-like quantities in our analysis, we consider 
each cluster from three different sightlines, namely along the x, y, 
and z axis of the simulation box. Unless noted, our analysis relies 
on the full cluster sample. 

A model galaxy is included in our samples if it lies within 
a threedimensional clustercentric distance of 3Mpc/ft. The clus- 
ter center is defined as the location of the central galaxy in the 
FoF halo, placed at the potential minimum of the FoF halo by 
the semi-analytic model. For the full and the SAM-V samples, 
we restrict model galaxies to a magnitude range of —19.0 < 
Mr < —15.2 mag. For the SAM-CP sample, the magnitude range 
is -19.0 <Mr < -16.7 mag. 

Some galaxies that lie beyond 3 Mpc/ft from the center may 
nevertheless be seen at a much lower projected clustercentric dis- 
tance. We test for such "interlopers" in the SAM-V clusters by ex- 
tracting all model galaxies with velocities in the range ±1800 km/s 
around the respective central cluster galaxyQ For the model clus- 



ters, we adopt a distance of 16.5 Mpc from a virtual observer ~ 
corresponding to the Virgo value - and select only galaxies with 
line-of-sight distances up to a factor of 1.5 larger or smaller The 
fainter or brighter absolute magnitudes that an observer would er- 
roneously derive when assuming them to be cluster members are 
taken into account for the selection. We expect that distinct struc- 
tures outside of this distance range, like the Virgo clouds mentioned 
in Section l2n could be identified by observers and would thus not 
be counted as cluster members in observational samples. 

At any projected clustercentric distance below 0.5 Mpc, the 
fraction of interlopers is 1% or less. For sightlines parallel to the 
x and y axis of the simulation box, the interloper fraction remains 
below 5% out to 1.5 Mpc, which is the completeness limit of the 
Virgo cluster sample. For sightlines parallel to the z-axis, this frac- 
tion remains below 10%. Between 1.5 and 2 Mpc, the interloper 
fraction does not exceed 10% for all sightlines. Therefore, our ne- 
glect of interlopers will not make a difference for the comparison 
of model clusters to the Virgo cluster For the Coma and Perseus 
clusters, which we consider up to 3 Mpc, the rising interloper frac- 
tion (reaching values of 30% and more) may seem relevant. On 
the other hand, we do coiTect the observed samples statistically 
for background galaxies, to which interlopers would be counted. 
Notwithstanding the uncertainty on the adopted background level 
(see Paper I) and possible effects of cosmic variance, the fact that 
the correction is done separately for red and for all galaxies re- 
duces a potential effect of interlopers on the colour-distance rela- 
tion, which is subject of our analysis in Sect. 14.41 

Finally, we remark that the stellar bulge-to-total mass ratio 
of our model galaxies is 0.03 or less for 50% of objects and 
0.10 or less for 75% of objects. Only 10% of model galaxies are 
bulge-dominated, i.e. have a bulge-to-total ratio of 0.5 or larger 
This seems consistent with the observational result that the vast 
majorit}lf| of these galaxies do not possess structural components 
with surface brightness profiles steeper t han a Sersic index of 2 
dOraham & Worlevl2008tlJanz et al.ll2013h . 



4 ANALYSIS 

The presence of significant tidal forces is a natural characteristic 
of ACDM cluster and group environments, both due to the over- 
all potential and to close encounters of subhaloes. The interplay of 
mass and gravitation cannot be avoided by the member galaxies, 
and causes significant mass loss of their subhaloes. However, we 
need to ask whether the stellar and gaseous components of a galaxy 
are necessarily affected when its subhalo loses mass. The link be- 
tween dark matter and baryons is provided by dedicated simula- 
tions of individual galaxies moving within the gravitational poten- 
tial of a galaxy cluster and experiencing encount ers with other sub- 
haloes (Moore et al. 1998; Gnedin 2003b; Mastro pietro et al.l200^; 
ISmith et al.l2010() . While depending on the specific trajecto ries and 
also o n the local structure of the gravitational potential dGnedinl 
l2003ah . the strongest effects on the baryonic configuration appear 



^ http://gavo.mpa-garching.mpg.de/MyMiIlennium/ 

^ This corresponds to the last snapshot of the simulation. In Paper I we used 
the second-to-last snapshot, at a lookback time of 263 Myr and a redshift of 
z = 0.020, which is closer to the redshift of the Coma and Perseus clusters. 
* This range is larger than two times the line-of-sight velocity dispersion 
of each model cluster sightline, except for one sightline to one cluster (Pa- 
per I). 



^ Note that our term "dwarf galaxies" refers to a simple magnitude se- 
lection of —19.0 < Mr < —15.2 mag, not to any selection in surface 
brightness. Observational samples of early-type galaxies in this range do in- 
clude a small fraction of objects classified as "low-luminosity elliptical/SO" 
galaxies (Janz & Lisker 2008), which differ from the more diffuse "dwarf 
elliptical/SO" galaxies. Having a small number of bulge-dominated galaxies 
in the sample thus seems reasonable from an observational point of view. 
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in cases with tiie strongest dark matter mass loss. We can thus as- 
sume that, statistically, the more dark matter was lost by a given 
subhalo over time, the more its stellar structure was tidally heated, 
leading to thickening or eve n destruction of disks and c ausing mild 
to str ong stellar mass loss dMastropietro et all 120051 : ISmith et^ 
I2OIOI) . In the course of our analysis, we will therefore focus on the 
mass loss experienced by the subhaloes in the Millennium-II sim- 
ulation. However, since mass loss itself does not directly indicate 
in which environment a galaxy resided, we additionally consider 
the total time that it has spent in massive haloes. This allows us to 
examine whether subhaloes that have suffered stronger mass loss 
have necessarily resided longer in high-mass environments. 



4.1 Infall time and mass loss 

In the upper panel of Fig. [T] we show the distribution of lookback 
times to when the model galaxies became a satellite for the first 
time, when they entered a halo with mass M > 10^^Mo//i, ahalo 
with mass M ^ 10^'^ Mq /h, and when they entered the cluster 
halo of which they are a member today. Since the simulation out- 
put is provided in discrete snapshots, we work with the simulation 
snapshot immediately before a galaxy appeared as member of the 
respective FoF halo for the first time. 

Almost 1600 galaxies (17.2%) are not yet a member of the 
cluster - in the sense of belonging to its FoF halo - but are located 
within 3Mpc/fe of the cluster center. These objects are assigned 
a lookback time of zero. They are no t included in the simila r dis- 
tributions of infall time presented by lOe Lucia et al.l ( l2012bl their 
Fig. 7), which are otherwise consistent with our distributions|f| 

While the distributions shown in Fig. [T] are a combination of 
15 clusters, a significant scatter between individual clusters exists, 
depending on which lookback time is chosen. The average look- 
back time to when 50% of today's dwarf galaxies had entered their 
cluster halo is 5.55 Gyr, with a rather large standard deviation of 
1.51 Gyr. This is mainly due to major merger events of the clus- 
ter haloes themselves, during which large groups or clusters are 
being accreted. These events thus set the cluster infall time for a 
significant fraction of today's galaxies. If we focus instead on the 
infall time into a halo with mass M ^ IO^^Mq/Zi - no matter 
whether or not it was the progenitor of today's cluster halo - the 
scatter is much smaller: the average lookback time is 7.33 Gyr with 
a standard deviation of only 0.55 Gyr between the clusters. For the 
infall time into a halo with mass M ^ 10^"^ Mq /h, the value is 
8.72 ± 0.45 Gyr. 

When choosing as reference the subhalo mass immediately 
before it became a satellite for the first time, a relative mass loss can 
be defined and compared for the different lookback times (lower 
panel of Fig. [Hi. Objects that have never been a satellite before the 
respective point of time are assigned a mass loss value of zero. This 
is why the peak in the figure at a value of zero increases with larger 
lookback time. Subhaloes that grew in mass since they first became 
a satellite appear with negative mass loss values in the figure, as 
indicated there. 

For a significant fraction of objects, a major part of the sub- 
halo mass had been lost already before entering today's cluster halo 
and even before entering a halo with 10^'^ Mq/H. Only when look- 
ing back to the infall into a halo with 10^"^ Mq / h, the fraction of 



^ Note that lGuo et all )201 lb and lOe Lucia et alj )2012bh use the most re- 
cent time when a galaxy has become a satellite, whereas we use the first 
time when it became a satellite. 
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Figure 1. Top: For the dwarf galaxies in the 15 most massive clusters of the 
Millennium-II simulation, we compare the lookback time to when a subhalo 
first became a satellite, when it entered a halo with mass M ^ 10^^ Mq / h, 
a halo with mass M 10^^ Mq/H, and today's cluster halo. Bottom: 
Here we compare the subhalo mass loss, relative to when it first became 
a satellite, for four different times: today (z = 0), when it first entered a 
halo with mass M 3= IO^^Mq/H, a halo with mass M ^ lO^^Me/h, 
and today's cluster halo. If one of these times is the same as the first time 
of becoming a satellite, the mass loss is set to zero. Negative mass loss 
values mean that the mass has increased since becoming a satellite. We 
assign orphan galaxies a subhalo mass corresponding to 19 particles, i.e. 
one particle below the resolution limit. 



objects that had already lost the majority of their mass before this 
event drops below 10% (black histogram in Fig. [T] lower panel). 
At that point of time, most subhaloes either experienced their first 
time of becoming a satellite - thus having a mass loss value of zero 
- or had even grown in mass since becoming a satellite. The latter is 
possible because many objects became first-time satellites at very 
early epochs (see the upper panel of the figure), when they were still 
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Figure 2. For the dwarf galaxies (white dots) in the 15 most massive clus- 
ters of the Millennium-II simulation, we show on the abscissa the look- 
back time to when a subhalo first entered a halo equal to or more massive 
than IO^'^Mq/Zi (right panels) and the subhalo mass loss since that time 
(left panels). On the ordinate, we show from top to bottom the time spent 
as satellite, the time spent in haloes with mass M ^ W^^Mq/H, with 
M ^ 10^^ Mq/H, and the time spent in today's cluster halo. Black cir- 
cles and error bars denote the median and the ±40% range in eight evenly 
spaced intervals. The black number in the top of each panel is the Spear- 
man's rank correlation coefficient multiplied by 100, calculated within the 
panel limits. 




1 2 
Projected clustercentric distonce / Mpc 

Figure 3. For the dwarf galaxies in the 15 most massive clusters of the 
Millennium-II simulation, we show the projected clustercentiic distance 
versus the logarithm of the local projec ted galaxy number density. The lat- 
ter is derived following IPresslerl l l 1 980h . by using the minimum radius of a 
circle that encompasses - in projection - the tenth nearest neighbour galaxy, 
counting all galaxies with Mr < —15.2 mag. The density is then calculated 
as 1 1 divided by the area of that circle in Mpc-^ . 



They obviously follow clear correlations, but also show signif- 
icant scatter. For example, the time spent in haloes with mass 
can be much shorter than the lookback time 
to the first infall into such a halo (right panel in second row from 
top), indicating that galaxies can pass through and escape a massive 
halo. Those (few) subhaloes that have grown in mass since their in- 
fall intoalO^^Mo/Zihalo can have spent very different amounts of 
time in such haloes, from almost zero to more than 7 Gyr (left panel 
in second row from top). Moreover, we emphasize that the scatter 
of the time spent in today's cluster halo is large (bottom panels): 
for many galaxies, this time just means the most recent stage of 
their evolution, not necessarily being representative of what they 
experienced before. 



very small in mass, and shortly afterwards they became again "cen- 
trals", i.e. the most massive galaxies of their respective haloes. This 
allowed them to continue growing significantly. Choosing as refer- 
ence the first time of becoming a satellite thus seems only moder- 
ately useful. Instead, our reference from now on will be the time 
when first enterin g a halo with mass M ^ 10^^ M^/h (also see 
iMcGee et al]|2009h . Only 2.2% of subhaloes possess a larger mass 
today than at this point of time, while the fraction is 17.3% when 
comparing to the mass at becoming a satellite. 

Since subhalo mass loss is caused by the tidal forces acting on 
a galaxy, we expect it to be a more direct proxy for the environ- 
mental influence - also on the baryonic part - than the lookback 
time to infall into a halo of a certain mass. Furthermore, the total 
time spent in massive haloes should also be somewhat closer re- 
lated to environmental influence than the lookback time to the infall 
event. In contrast to the mass loss, which depends on the orbit and 
the occurrence of encounters, the time spent in massive haloes can 
tell whether the progenitors of different galaxy populations have 
evolved in different environments, thus providing complementary 
information. 

The three quantities are compared to each other in Fig. |2] 



4.2 Present-day location versus environmental influence 

The quantities of the subhalo population that are most directly ac- 
cessible to an observer are the projected position and line-of-sight 
motion. We will address the line-of-sight velocity in Section [43] 
but focus here on location, i.e. projected distance from the respec- 
tive cluster center In addition, we also consider the commonly 
used quantity local density, i.e. the projected number density of 

i ;alaxies calcula ted from the area encompassing the 10th neighbour 
DressleJll980l) . Both are compared in Fig.|3] showing a clear cor- 
relation, but with significant scatter Part of the scatter is due to 
cluster substructure: a subclump with relatively high local density 
may be located comparatively far away from the cluster center. This 
is the case, for example, with the Virgo cluster and its M 49 subclus- 
ter, which is more than 1 Mpc away in projection from the central 
Virgo galaxy M 87. 

In Fig.|4]we show how both clustercentric distance and local 
density correlate with subhalo history, i.e. with lookback time to 
infall, time spent in haloes with M ^ 10^^ Mq / h, and mass loss. 
All correlations are similarly strong for clustercentric distance as 
for the logarithm of local density, judging from the visual impres- 
sion as well as from the correlation coefficient given in each figure 
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Figure 4. For the dwarf galaxies (white dots) in the 15 most massive clus- 
ters of the Millennium-II simulation, we show on the abscissa the projected 
clustercentric distance (left panels) and the logarithm of the local projected 
galaxy number density (right panels). On the ordinate, we show from top to 
bottom the lookback time to when a subhalo first entered a halo with mass 
M ^ IO^^Mq/Zi, the lookback time to when it entered today's cluster 
halo, the time spent in haloes with mass M > IO^^Mq/H, the subhalo 
mass loss since it first entered a halo with mass M ^ 10^-^ Mq/H, and the 
mass loss since it entered today's cluster halo. Black circles and error bars 
denote the median and the ±40% range in eight evenly spaced intervals. 
The black number in the top of each panel is the Spearman's rank correla- 
tion coefficient multiplied by 100. calculated within the panel limits. 



panel. Objects located at small clustercentric distances, and like- 
wise in regions of high local density, experienced on average early 
infall, have spent a long time in massive haloes, and have suffered 
strong mass loss. However, the scatter is large, clearly illustrating 
that the past environmental influence can not simply be read off 
from the (projected!) clustercentric distance today. There is more 
than 50% overlap between the subhaloes at less than 0.5 Mpc from 
the cluster center and those between 1 and 2 Mpc in terms of the 
time spent in haloes with AI 10^'^ AIq / h, the mass loss expe- 
rienced in such haloes, and also the lookback time to infall into 
today's cluster. Only part of this overlap is caused by projection 
effects — even when considering threedimensional clustercentric 
distance, the statement still holds. We conclude that the subhalo 
populations in the cluster center and the outer cluster regions do 



show systematic differences in their histories on average, but that 
two subhaloes located at the same clustercentric distance or local 
density today may have experienced strongly different environmen- 
tal influence in the past. Two subhaloes located at very different 
clustercen tric distances o r local densities may ha ve had similar his- 
tories (cf. lGnedinir2003al : lDe Lucia et al. I l2012bh . In the following 
section, we will attempt to map these subhalo properties to the ob- 
served galaxy population. 

To ensure that our selection of model galaxies in r-band mag- 
nitude instead of stellar mass does not introduce a bias, we examine 
the above relations with clustercentric distance for model galax- 
ies in the following stellar mass intervals: 0.1 - 0.4 x 10** Mq, 
0.4 - 2.0 X 10® Mq, 2 - 10 X 10* Mq, 10^ - lO^^M©, and 
10^° - 10^^ Mq. We find that the radial trends, as well as the 
scatter of values, are remarkably similar in all but the most mas- 
sive interval. For the lookback time to entering a halo with M ^ 
10^^ M0//1, the mass loss since then, and the time spent in such 
haloes, the median values differ by less than 10% at projected clus- 
tercentric distances below 1 Mpc and less than 20% out to 3 Mpc. 
The most massive interval differs by less than 20% below 1 Mpc 
and up to 70% out to 3 Mpc. The ±40% ranges differ by less than 
20% except for the mass loss below 1 Mpc, where the range can be 
different by a factor of two. In the most massive interval outside of 
1 Mpc, the ±40% range of mass loss already includes objects that 
have grown in mass (i.e. with negative mass loss values), thereby 
increasing the range significantly. However, for low-mass galax- 
ies, this test confirms that our analysis is not biased by our specific 
magnitude selection of model galaxies, which is tied to the obser- 
vational samples. 



4.3 Morphology-distance relation 

We now assume that, for dwarf galaxies, stronger environmental in- 
fluence in the past led to early-type morphology today. The present- 
day late-type galaxies would therefore be those on which environ- 
ment had the weakest effect. We thus attempt to assign subhalo 
populations with different mass loss histories to the observed pop- 
ulations of galaxies with different moiphology. The latter are taken 
from the Virgo cluster and are subdivided into three classes: (i) el- 
liptical and dwarf elliptical galaxies without disk features or blue 
central regions; (ii) lenticular galaxies, Sa-type spirals, dwarf ellip- 
ticals that exhibit disk features and blue central regions, and tran- 
sition types between dwarf elliptical and irregular galaxies; (iii) all 
remaining spiral galaxies, as well as irregular and blue-compact 
dwarf galaxies. When considering only Virgo galaxies up to a pro- 
jected distance of 1.5 Mpc from M 87 (which we assume to be the 
cluster center), these classes comprise 50.1%, 21.9%, and 27.9% 
of galaxies, respectively. Their distribution with clustercentric dis- 
tance is shown in the top pan el of Fig. 5] refl ecting the well-known 
morphology-density relation ( lDressleill980l) . 

The above percentages are now used to subdivide the model 
galaxies by their subhalo mass loss, and alternatively by the time 
spent in haloes with AI > 10^^Mq//i, assuming these quanti- 
ties represent the strength of environmental influence. The resulting 
distributions with clustercentric distance are shown in the two mid- 
dle panels of Fig. |5] They look similar to the observed relation of 
galaxy morphologies, but have a larger contrast between the inner 
and outer populations. Subhaloes that experienced strong mass loss 
or spent a long time in massive haloes dominate in the center, but 
already beyond 0.75 Mpc they are outnumbered by those with weak 
mass loss or a short time spent in massive haloes. In contrast, the 
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Figure 5. Clustercentric distribution of galaxies and subhaloes. The top 
panel shows the Virgo cluster moiphological types. These type percentages 
- using only galaxies up to 1.5 Mpc projected distance - are then used to 
subdivide the subhalo population accordingly. The second and third panel 
from top show the resulting distributions when subdividing by the time 
spent in massive haloes, or by the mass loss. For the lower panel, we as- 
sume the opposite mapping of observed and model galaxies, to illustrate 
that this would lead to a discrepant distribution. We use the SAM-V sample 
of model clusters for this figure (see Section|3]. 



elliptical galaxies in the Virgo cluster dominate over the late types 
out to 1.5 Mpc. 

The above assignment of percentages assumes that strong 
mass loss correlates with early-type morphology. When we as- 
sume instead the opposite, namely that late-type galaxies are those 
that experienced the strongest mass loss, the lower panel of Fig. |5] 
shows the resulting relation with clustercentric distance. This sce- 
nario would lead to a relation that is opposite to the observed 
morphology-distance relation, with a strongly increasing elliptical 
fraction towards the cluster outskirts. This illustrates that our above 
assumption was reasonable, and that one cannot simply assume any 
arbitrary correlation between dark matter mass loss and baryonic 
morphology. If we assumed that no correlation exists, this would 
lead to a scenario in which any given galaxy type would have to 
cover the full range of mass loss values. For example, the popu- 
lation of late-type, star-forming galaxies would have to be com- 
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Figure 6. Similar to Fig.|5] but now using the morphological percentages 
of each individual radial bin of the Virgo cluster distribution to subdivide 
the model galaxy population of that bin by mass loss (top panel). Points and 
eiTor bars denote the median and the ±40% range. The bottom panel shows 
the resulting distribution of the time spent in haloes with ^ IQi^'^ Mq / h. 
We use the SAM-V sample of model clusters for this figure (see Section[3j. 



posed partly of galaxies whose subhaloes have lost almost all of 
their mass, but also of galaxies whose subhaloes have experienced 
almost no mass loss. All of these would intriguingly have to show a 
very similar appearance today, since all belong to the same galaxy 
type. This appears rather unlikely, given that N-body simulations of 
tidal forces in galaxy clusters show that objects losing a major frac- 
tion of their dark matter are affected in their st ellar configuration 
as well ( lGnedirj|2003al : lMastropietro et al.l2(j05l) . 

When we subdivide the model galaxies of the SAM-V sam- 
ple by their subhalo mass loss as outlined above, the median 
g — r colour of the early-type analogues turns out to be 0.70, which 
is slightly redder than the intermediate-type analogues (0.66) and 
clearly redder than the late-type analogues (0.55). This is a conse- 
quence of the implementation of environmental effects in the semi- 
analytic model (see Paper I an d iGuo et al.l20T l| for details): objects 
with stronger mass loss have become satellites of massive haloes 
at earlier times (Fig. |2j. While it lends support to our approach, it 
needs to be remarked that model colours are too red as compared to 
observations (see the analysis and discussion in Paper I). The Virgo 
cluster early-type dwarfs have a median g — r colour of 0.60, the 
intermediate types have 0.60 as well, and the late types have 0.39. 

Instead of using only one set of population percentages, de- 
rived from all galaxies up to 1.5 Mpc from the center, we can map 
model galaxies to observed ones for each clustercentric distance 
bin individually. This is shown in the top panel of Fig. [6] using 
mass loss to subdivide the populations. As a result, we can see that 
there exists a relation of mass loss and distance within each class 
of galaxies. This is a direct consequence of the fact that the rela- 
tion with distance in Fig.|5]is steeper for model galaxies than for 
observed ones when using & fixed value to separate galaxy popula- 
tions: the decline within the populations was not taken into account. 
It is most pronounced for the late-type analogues, which reach a 
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Figure 7. Red dwaif galaxy fractions (black line) and subhalo fractions with 
strong mass loss (grey line and shaded area) are shown with respect to pro- 
jected clustercentric distance for observed clusters (Virgo, Perseus, Coma) 
and their corresponding model cluster samples (SAM-V for Virgo, SAM- 
CP for Coma and Perseus, see Section[3] . Red galaxies are defined accord- 
ing to eq. [T] following Paper I. The shaded areas indicate the minimum- 
to-maximum range of the different model clusters and their different pro- 
jections. Subhaloes with strong mass loss are selected similar to Fig.|5] we 
apply the percentage of red galaxies inside a projected distance of 1.5 Mpc 
(Virgo) and 2.0 Mpc (Perseus and Coma) to the model samples. The Virgo 
sample becomes incomplete beyond 1 .5 Mpc, indicated by the vertical dot- 
ted line. The Perseus sample is statistically coiTected for incomplete cov- 
erage (Sect. 12.3) , but this becomes unreliable beyond 2.5 Mpc due to too 
small coverage. 



median mass loss value of 50% in the inner cluster regions, but 
drop to less than 15% in the outskirts. This could indicate a mass 
loss threshold for the onset of noticeable effects on the stellar con- 
figuration. Without a major loss (> 50%) of subhalo mass, galaxies 
were able to keep their late-type appearance. 

In addition to mass loss, we show in the bottom panel of Fig.|6] 
the time spent in massive haloes that results from the above sub- 
division by mass loss. Again, there is a significant trend within 
each class, especially for intermediate- and late-type galaxies. 
Intermediate-type galaxies in the innermost radial bin have spent 
a similar time in massive haloes as elliptical galaxies located in the 
outskirts. The same is true for late-type galaxies in the innermost 
radial bin as compared to intermediate-type ones in the outskirts. 

However, we have so far considered the combination of all 
model clusters of the SAM-V sample, and all their projections. If 
the Virgo cluster is not a typical cluster, the actual distribution of 
mass loss with clustercentric distance may be different than shown 
in Fig.[6] We have inspected the individual projections of the twelve 
SAM-V clusters, and found that in a few of the 36 cases, the relation 
within a given subclass almost disappears. Moreover, the stellar and 



gas structure of a galaxy may be governed by further parameters 
in addition to subhalo mass loss, which could help understanding 
why late-type galaxies in the center have mass loss values similar 
to intermediate-type galaxies in the outskirts. These aspects illus- 
trate how complex it is to interpret the observed galaxy populations 
correctly. 

4.4 Colour-distance relation 

A similarly complete galaxy catalog as for the Virgo cluster does 
not exist for the C oma and Perseus clusters (only for the Coma cen- 
ter, see iMichard & AndreonllZOOsI) . but an analogous comparison 
between observed and model galaxies can be done based on galaxy 
colour, as shown in Fig. [7] While colour is not a direct proxy for 
morphology, galaxies in most clusters show a colour-distance rela- 
tion similar to the morphology-distance relation. This can be seen 
in the top and bottom panels of the figure: the red galaxy fraction 
of the Virgo and Coma cluster declines with increasing distance 
(black line). Red galaxies are selected by eq. [T] after applying a 
k-correction to Coma and Perseus galaxies (Sects. l2!2l and l2.3T ). In- 
terestingly, the red galaxy fraction of the Perseus cluster does not 
decrease at all when going outwards. This shows that significant 
scatter exists among the properties of present-day massive galaxy 
clusters (also see Paper I), and an individual cluster can not neces- 
sarily be taken as representative of the majority of clusters. 

To obtain distributions of subhaloes with strong mass loss, we 
apply for each observed cluster the overall percentage of red galax- 
ies to its corresponding model cluster sample, analogous to Fig. [5] 
The resulting subhalo distributions are shown with grey lines, and 
the shaded area indicates the full range of distributions that occur 
among the various projections of the different model clusters. For 
the Virgo cluster, a similar trend is seen as for morphology: the 
fraction of subhaloes with strong mass loss decreases more strongly 
than the observed fraction of red galaxies. Nevertheless, the Virgo 
cluster still lies within the range of the different model cluster pro- 
jections, same as the Coma cluster in most radial bins. The Perseus 
cluster with its flat observed distribution falls outside of the model 
cluster range in half of the radial bins. However, it needs to be em- 
phasized that the SAM-CP sample consists only of the three most 
massive model clusters. With a larger model sample - requiring a 
cosmologial simulation with larger box size at the Millennium-II 
resolution - the range of distributions may become larger. 

4.5 Mass loss history of galaxy populations 

To what extent is the present-day population of late-type galaxies 
in clusters comparable to the former progenitors of today's early 
types? We address this question by comparing the mass loss dis- 
tributions of these populations with each other, assuming that sub- 
halo mass loss serves as a proxy of the environmental influence on 
the (baryonic) galaxy. Following the subdivision of Virgo cluster 
galaxy types described in Section [43] we select subhaloes repre- 
senting elliptical, intermediate, and late-type galaxies in the model 
clusters, and show their mass loss distributions in Fig.|8] In the top 
panel of the figure, we subdivide the three galaxy types based on 
the time spent in haloes with ^ W^^ Mq/H, while in the bottom 
panel we subdivide them based on the mass loss itself. As before, 
mass loss is calculated relative to the subhalo mass immediately 
before entering a halo with M ^ W^^ Mq / h for the first time. 

The mass loss distribution of the present-day late-type galax- 
ies is illustrated by the grey-shaded histogram in the figure. A sig- 
nificant number of objects - note the logarithmic ordinate - have 
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Figure 8. Shown are the distributions of mass loss since entering a halo 
with ^ IO^^Mq/Zi for the model galaxy populations defined as in Fig. [5] 
As in that figure, the model populations ('elhpticar, 'intermediate', and 
'late-type') can be defined either based on the time spent in massive haloes 
(top panel) or on the mass loss they experienced (bottom panel). For the 
'elliptical' model galaxies only, we show how the mass loss distribution 
evolves with redshift (coloured lines), to find out at which redshift (if any) 
it has been compai'able to the present-day late types. The total numbers 
of intemiediate- and of late-type galaxies were scaled to match that of the 
elliptical galaxies. 



been assigned zero mass loss, since they are not FoF-member of 
their respective clustei[3| and have never been a FoF-member of a 
halo with M ^ 10^^ Mq/Zi. Negative mass loss values mean that 
the mass increased since infall. In the top panel of the figure, some 
late-type galaxies even experienced a complete mass loss, i.e. a 
value close to 1. In the bottom panel, this is not possible by con- 
struction: mass loss serves as the basis for the subdivision itself, 
therefore late-type galaxies reach only to values of about 0.6. Note 
that the numbers of intermediate- and late-type galaxies were nor- 
malised to that of the elliptical galaxies. 

For the present-day elliptical galaxies (black histogram), we 
show the redshift evolution of the mass loss distribution with the 
coloured histograms. With increasing redshift, the number of galax- 



^ We remind that we include in our analysis all subhaloes located within 
3 Mpc//i of the cluster center, no matter whether or not they are counted 
as members of the FoF-halo. 
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Figure 9. For the dwarf galaxies (white dots) in the 15 most massive clus- 
ters of the Millennium-II simulation and located within 0.5 Mpc of the clus- 
ter center, we show on the abscissa the absolute line-of-sight velocity rel- 
ative to the cluster center (left panels, combining all three projections) and 
the absolute threedimensional velocity relative to the cluster center (right 
panels). On the ordinate, we show from top to bottom the lookback time 
to when a subhalo first entered a halo with M ^ 10^^Mq//i, the look- 
back time to when it entered today's cluster halo, the time spent in haloes 
with M ^ IO^^Mq/Zi, the subhalo mass loss since it first entered a halo 
with M > IO^^Mq/Zx, and the mass loss since it entered today's cluster 
halo. Black circles and error bars denote the median and the ±40% range 
in eight evenly spaced intervals. The black number in the top of each panel 
is the Spearman's rank congelation coefficient multiplied by 100, calculated 
within the panel limits. The restriction to a clustercentric distance of less 
than 0.5 Mpc is applied to the projected distance for the left panels and to 
the threedimensional distance for the right panels. 



ies that had already experienced strong mass loss is lowered. At 
the same time, the peak at the value of zero grows higher, in- 
dicating that these objects were not yet member of a halo with 
M ^ IO^^A/q/Zi at the respective redshift. However, the distri- 
bution never looks the same as that of the present-day late types. 
These have a distribution that decreases at large mass loss values 
(upper panel of Fig. [8), while the distribution of the ellipticals' pro- 
genitors either increases or is nearly flat. If one assumes that only 
strong subhalo mass loss - e.g. more than 50% - had a noticeable 
effect on the stellar configuration of the galaxies, then only at red- 
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shifts between 2 and 3, the elhpticals' progenitors become similar 
to the late types of today. 



4.6 Dependence on velocity 

In addition to the clear correlations of projected position with infall 
time and mass loss, we investigate whether these quantities show 
an additional correlation with velocity. For this purpose we con- 
sider the absolute value of the threedimensional velocity relative 
to the central cluster galaxy, as well as the absolute value of the 
respective line-of-sight velocities along each of the three axes of 
the simulation box. The line-of-sight velocities do not show any 
significant correlation with the lookback time to when a subhalo 
first entered a halo with A/ ^ 10^'^ Mq / h, the lookback time to 
when it entered today's cluster halo, the time spent in haloes with 
M ^ IQ^'^Mq/H, the subhalo mass loss since it first entered a 
halo with M ^ IQ^^ Mq /h, nor with the mass loss since it entered 
today's cluster halo. 

However, when restricting our analysis to galaxies located 
within a projected clustercentric distance of 0.5 Mpc, i.e. in the 
cluster center, weak correlations become apparent (Fig. |9] left pan- 
els). While the trends are smaller than the scatter of values, and 
the individual correlation coefficients are small, the relations give 
a consistent picture: galaxies moving with larger relative veloci- 
ties had a somewhat later infall, have spent less time in massive 
haloes, and have experienced less mass loss. In fact, these corre- 
lations are smeared out by projection effects — they are clearly 
stronger when using threedimensional quantities (right panels of 
Fig.|9). This probably reflects the fact that galaxies that entered the 
gravitational potential of the halo at a later time, when the cluster 
had already grown to a larger mass, acquired larger velocities until 
reaching the center. 

The comparison between projected and threedimensional 
quantities also serves as example for the effect of being in the ob- 
server's situation and having only projected quantities available: 
strong correlations may be somewhat weakened, and weak corre- 
lations may become insignificant. We point out that no further or 
stronger correlations become apparent when using not the absolute 
values of lookback time and mass loss, but their residuals about the 
relation with projected distance. 

Our analysis of subhalo velocities cannot fully explain the ob- 
served strong correlation between line-of-sight velocity and shape 
of nucleated dwarf elliptical galaxies in the Virgo cluster center 
jLisker et al] l2009t) . The finding that those galaxies with lower 
line-of-sight velocities have significantly rounder shapes had been 
interpreted such that these have resided in the cluster since a 
longer time, are therefore on more circularised orbits (see also 
iBiviano & Poggiantill2009l) and have suffered more environmental 
influence. However, while we do find weak correlations of line-of- 
sight velocity with subhalo mass loss and infall time, these are far 
from being as clear as the observed correlation with shape, which 
remains to be understood. 



4.7 Dependence on mass 

Now we address the question whether the differences between the 
galaxy populations in the inner and outer cluster regions become 
more or less pronounced with increasing cluster mass. Fig.llOlcom- 
pares the lookback time, the time spent in massive haloes, and the 
mass loss of the inner (white) and outer (black) model galaxy pop- 
ulation of each simulated cluster. The inner and outer red galaxy 



.1 § 

So 
n -- 

^ s 



i_ o 
c 2 



ul O 
£ All 



10 

5 


10 



l/l CM 



0.5 



p 0.5 



outer 
(Mpc) 




H 1 1 1 1 1 1 1- 



H 1 l-^f — I 1 1 1 1 1- 



Virgo Perseus Coma 



:|> o 
outer 

(Rvir) I 



H 1 1 1- — I 1 1 l-^l 1- 



H 1 1 1 1 1 1 1 1 1- 



Virgo Perseus Coma 



14.5 15 14.5 15 

log( Virlal Mass of Cluster / Mg ) 

Figure 10. Median and ±40% range of model quantities for the inner 
(white) and outer (black) subhalo populations of the individual model clus- 
ters, shown with respect to cluster mass. The bottom panel shows the red 
galaxy fractions of the observed clusters. In the left panels, "inner" and 
"outer" is defined in units of Mpc: less than 0.5 Mpc projected clustercentric 
distance, and the range 1.0 — 1.5 Mpc, respectively. In the right panels, we 
define "inner" and "outer" in units of virial radii: less than 0.33 -Rvir and 
the range 0.67 — 1.0 -Rvir, respectively. For the Virgo cluster, these defini- 
tions are the same, since the virial radius is assumed to be -Rvir = 1-5 Mpc 

i McLauehlin 1999). For the Perseus cluster we adopt -Rvir = 1.8 Mpc 

— (f i 

(estimated from extrapolating the mass profile of Eyles et al. 199 11, see Pa- 
per I) , and for the Coma cluster we use -Rvir = 2.8 Mpc iLokas & MamoiJ 



fractions of the observed clusters are shown for comparison. Note 
that the left panels of the figure distinguish between "inner" and 
"outer" regions in units of Mpc, choosing projected clustercentric 
distances less than 0.5 Mpc and the range 1.0 — 1.5 Mpc, respec- 
tively. In contrast, the right panels are based on units of virial 
radius 7?vir, with less than 0.33 -Rvir for "inner" and the range 
0.67 — 1.0 -Rvir for "outer", respectively. 

The lookback time to when a model galaxy first entered a halo 
with M ^ 10^^ Mq/H is always high for the inner population. A 
trend is seen for the outer population, when defined in units of Mpc, 
that lookback time is increasing with cluster mass (top left panel of 
Fig. not . However, this is mostly an effect of choosing absolute 
units: no significant trend is seen when defining the outer popula- 
tion in units of -Rvir (top right panel). The situation is the same for 
the time spent in massive haloes (second row from top), although 
a slight tendency may be present in units of -Rvir. When focusing 
on mass loss (third row from top), again the inner populations are 
all very similar and have a mass loss close to one. The outer popu- 
lations, in units of Mpc, show a clear trend of increasing mass loss 
with cluster mass, and a weak tendency in units of -Rvir. 
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The red fractions of observed galaxies show no clear trend 
for the inner populations. The outer populations follow a slight in- 
crease with cluster mass when defined in units of Mpc, and no trend 
when defined in units of 7?vir- This may appear similar to the model 
galaxies, but is of course subject to small number statistics. 

In our previous diagrams with clustercentric distance, we have 
combined the different model clusters and their projections in units 
of Mpc. Therefore, the differences between inner and outer popu- 
lations may have been smoothed out a bit, as compared to when we 
would have used units of J?vir. On the other hand, it is not trivial 
to estimate the virial radius of a cluster from observations, which is 
why we decided to present our diagrams in units of Mpc instead. 



5 DISCUSSION 

5.1 Limitations of linking sublialoes with their baryons 

Despite the exquisit mass resolution of the Millennium-II simu- 
lation, the number of particles that represent a low-mass subhalo 
is rather small: about 1500 particles for a subhalo with M = 
IO^^Mq//!.. Thus, the actual dynamical interaction of tidal field 
and subhalo, or of two subhaloes, cannot be properly described 
by a cosmological simulation itself. Moreover, many of the model 
galaxies in the central cluster regions are orphans, i.e. their sub- 
haloes dropped below the resolution limit of 20 particles and are 
only tracked by the semi-analytic model. Observations indicate, 
howeve r, that t hese objects still hold a significant amount of dark 
matter jPennv et al. 2009). Therefore, the mass loss, as well as the 
position and motion of low-mass subhaloes, are merely approxima- 
tions. Furthermore, dark matter may couple to baryonic processes 
that change the gravitational potential locally, like gas loss induced 
by supemovae or by ram pressure. These can affect the dark matter 
profile and distribution (Govemato et al. 2010; Smith et al. 2012i). 

Subhalo mass loss begins to occur earlier and already in a less 
strong tidal field than major changes of the baryonic configuration. 
This is due to the large extent of the initial halo, which then loses 
its o uter parts and gets truncated when entering a larger potential 
dOned in 2003b; Villalobos et al. 2012). Therefore, when inteipret- 
ing subhalo mass loss with regard to observable effects, a major 
fraction of the subhalo can be lost without affecting the baryonic 
galaxy. This is in line with the majority of late- type analogues hav- 
ing lost < 50% of their subhalo mass (Fig.|6}: not much harm has 
yet been done to their stellar disks. 

While our study focuses on the mass loss of dark matter sub- 
haloes caused by tidal forces, the strongest ram pressure occurs 
in the same environments: the highest hot ga s densities are found 
in m assive galaxy clusters jMohr et ai]|l999l : iHelsdon & Ponmai] 
12000) . To what extent ram pressure stripping contributes to quench- 
ing the star formation activity of galaxies falling into a group or 
cluster is still a matter of debate ('Goto"2005';'McCarthv et al."200_^; 
Book & Benson 2010; Weinmann et al. 2010; Roediger et al. 2011i; 
Bahe et ali2013l ; lBosch et alj2013l) . It is also noteworthy that there 



jLisker et~^ 20061: TuUv & TrenthamI 2008 ) and contain g as and 



exists a significan t scatter in the intragroup medium properties 
of ga laxy groups jHelsdon & Ponmai] |2000| ; lOsmond & Ponmai] 
12004) . Studies aiming to fully simulate the origin of today's dwarf 
galaxy population thus need to model the variety of environments 
that affected the galaxies over their lifetime (see Fig.[T}. 

Simulations and models of envirormiental processes would 
also need to account for spatially resolved galaxy proper- 
ties. While having early-type morphology, dwarf galaxies in 
groups and clusters may still form stars in their inner regions 



dust I di Serego Alighieri et al.ll2007l; lOe Looze et al.ll2010^ . Even 
among those not forming stars anymore, the majority of bright 
early-type dwar fs in the Virgo clust er exhibit a young central stel- 
lar population (Paudel et al.l 2010ah. The blue ultraviolet— optical 
central colours (Boselli et al.l l20o'8ah o f galaxies with overa ll red 
colours and old stellar population ages ('Roediger et al. 201 1) may 
indicate a recent ram-pressure stripping event ( Boselli et al. 2008al) 
or even the reaccretion of gas dHallenbeck et al.ll2012l) . 

These observations seem to indicate a recent arrival of many 
early-type dwarf galaxies to the cluster environment, which could 
mean a discrepancy with the ACDM prediction that they are old 
objects ( iBoselli et al] l2008al : cf. Figs. □ and [ej. However, it is 
not straightforward to translate stellar population ages into sub- 
halo histories . When using the ram pressure stripping criterion of 
iGunn & Gotll ( 1 19721) as ap proximation, i t can be shown that the 
Virgo intracluster medium l lVollmedl2009[) would not be able to re- 
move gas from the centers of most bright early-type dwarf galaxies 
(Mr < — 17mag)[f| Only for fainter galaxies should its ram pres- 
sure seriously affect all parts of the galaxy, thus predicting that the 
Next Generation Virgo Cluster Survey (Ferrarese et al. 2012) will 
not find faint early-type dwarfs with blue central regions. 

The above considerations emphasize the complexity of the 
real situation, involving gas, stars, their time evolution, and their 
dependence on local conditions. The subhalo distributions and 
histories that we analysed constitute the underlying cosmological 
framework — not more and not less. 

A further caveat lies in the applicability of current semi- 
analytic models of galaxy formation in cosmological volumes, 
which are still facing difficulties in reproducing the properties of 
low mass galaxies. For example , they predict too little late evo- 
lution in their number density JWeinmann et al.l |2012|) . too red 
colours in galaxy groups (Paper I), as well as stellar metallicity and 
age distributions that do not agree with observations dPasguali et al.l 
2010). These aspects are one reason why we primarily rely on dark 
matter subhalo distributions for our analysis, instead of using the 
"observables" provided by the semi-analytic model for our com- 
parison with observations. 



5.2 Correlations with clustercentric distance 

The anticorrelations of mass loss and lookback time to infall with 
clustercentric distance (e.g. Fig. l4l lSmith et al.l2012dh seem to pro- 
vide a n atural explanation fo r the observed morphology-distance 
relation jBinggeli et al.|[l987l) . Spending a longer time inside mas- 
sive haloes means to experience stronger tidal forces and presum- 
ably also stronger interaction with the intracluster medium, leading 
to an earlier and more efficient transformation of the m orphological 
properties an d quenching of the star formation activity jSmith et al.l 
l2008ll2012c^ . 

Early-type dwarf galaxies al so show a morphology-distance 
relation of their own subtypes jLisker et al.. i2007i) , as well as 



* At a clustercentric distance of 0.5 Mpc , the de nsity of the Virgo in- 
tracluster medium as modeled by IVoUmed | |2009|) provides a ram pres- 
sure of 10^^^'*Nm~^. We select Virgo early-type dwarf galaxies with 
Mr < —17 mag and projected axis ratio above 0.85, and use ugriz- 
photometry to appr oximate their stellar surface mass density (based on 
iHansson et alj2012 ) at one exponential scale length. When assuming a gas 
surface density of one tenth of the stellai' surface density and a relative ve- 
locity of 1000 km/s, we obtain estimates for the restoring force per unit area 
between IQ-i^ ^ and IQ-l'^ O Nm'^. 



© 2013 RAS, MNRAS QOQ.fTHTTl 



History of cluster dwarf galaxy populations 13 



a dependency of t he colour-magnitude relation on local density 
jLisker et alj 12008*). Those early-type dwarfs with signatures of 
disks and/or without a bright stellar nucleus are preferentially found 
outside of the cluster core. These also exhibit younger ages of 
their stella r populations, deduc ed from stronger Balmer absorp- 
tion lines jPaudel et al.l |20103 for the Virgo cluster), fro m nar- 
rowband Stromgren photometry dRakos &, Schomber t 2004 for the 
Coma and Forna x clusters), and from u ltraviolet— optical c olours 
jKim et alj2O10l for Virgo). Accordinglv. lSmith et al.l ( l2008h found 
that red-sequence dwarf galaxies in the Coma cluster outskirts 
show stronger Balmer lines than those in the center. This can 
be understood in the light of our analysis: even within a given 
subpopulation, a radial trend exists for mass loss and for the 
time spent in massive ha loej^ (Fig. |6l cf. ISmith et^l2006l and 
Ivon der Linden et"^l2010i for bright galaxies). 

On the other hand, the various structural subgroups of Virgo 
early-type dwarf s that were identi fied in the comprehensive near- 
infrared study of Ijanz et al. l ll20Tl do not exhibit clear trends with 
respect to clustercentric distance. IRvs et all ( I20T2I) point out that 
the individual histories of the galaxies should lead to a significant 
spread in their characteristics at a given distance from the clus- 
ter center. This is in line with the substantial scatter in subhalo 
mass loss and time spent in massive haloes when considering the 
combined population of dwarf elliptical analogues and intermedi- 
ate types (Fig.[6l(. 

After it had become clear that a major fraction of early- 
type dwarf galaxies exhibit significant r otational velocities 
dSimien & Prugnielll2002l; IChilingarianll2009l) , | Toloba et al. Il l200^ 
reported a tendency for them to show increasing rotational support 
with increasing distance from the Virgo cluster center While this 
would seem complementary to the above trends of morphology and 
stellar population s, the kinematic al diversity of early-type dwarf 
galaxies is large ( IRvs et alj |2 01 2h and published samples are still 
incomplete. Again, a significant scatter is expected from our anal- 
ysis, since clustercentric distance cannot be directly translated into 
the specific history of a subhalo. 

Those early-type dwarf galaxies that reside in cluster cores 
have experienced high-density environments already at early 
epochs (Fig.[8l(. For them, more intense early star formation could 
have caused the larger fraction of s t ellar mass in the form of 
globular clusters jMoore et all bOOd IPeng et al.l l2008h as com- 

f iared to galaxies in today's clus ter outskirts or even the field 
Sanchez- Janssen & Aguerrill2012l) . Interestingly, the fraction of 
early-type dwarfs with bri ght s tellar nuclei that exceed typical glob- 
ular cluster lumin osities t Cote et al. 200 6j) is also much larger i n 
the cluster center d Ferguson & Sandagall989l iLisker et al.ll2007h . 
We can thus speculate that the efficient early formati on of globular 
clusters (and possibly ultra-compact dwarf galaxies, iMieske et al.l 
l2012i) in those regions was paralleled by an efficient formation of 
such stellar nuclei. This shows that, on top of the subhalo distri- 
bution and evolution, it is necessary to take into account the envi- 
ronmental and time dependence of baryonic processes to obtain a 
consistent picture of the formation of today's dwarf galaxies. 



5.3 The dwarf galaxy population in a ACDM universe 

Fro m their analysis of the dwarf ga laxy population in nearby clus- 
ters, ISanchez- Janssen et ai] (l2008f) concluded that the red dwarfs 
in the central cluster regions could be related to the population of 
blue dwarf galaxies observed in high-redshift clusters. However, 
at the brighter dwarf magnitudes, the red sequence was even re- 
ported to be es tablished and well pop ulated in clusters up to red- 
shift z ^ 1.3 l lAndreonlliooel 120081) . Many of the red-sequence 
galaxies in cluster cores experienced environmental effects already 
at hig h redshifts (Figs. [7] & [8) and entered the cluster halo very 
early dOe Lucia et al.ll2012bl : lsmith et al.ll2012dl) . 

On the other hand, clusters grow by accretion of field and 
group galaxies (cf. lAdami et al. ^ '2005*). Out of the 7673 model 
dwarf galaxies that are FoF members of the 15 most massive 
Millennium-II clusters at redshift zero, two thirds (66.9%) joined 
their present-day cluster halo at redshifts z < 1, and sti ll more than 
one th ird (36.4%) joined at redshifts z < 0.5 (see also lSmith et al.l 
'2009'). The majority of galaxies (57.2% ) were accrete d as satellites 
of a gr oup or cluster (also see iMcGee et al. 2009. De Lucia et al.l 
l2012bl . and lSmithetai]|2012dh . Envirormiental effects on galax- 
ies in groups before cluster accretion - so-called pre-processing - 
are therefore expe c ted to be of high relevan c e, as was noted also 
by iDe Lucia et all j2012bl) . Ivillalobos etai] i l2012h showed from 
simulations that stellar disks can indeed be significantly affected 
(thickened, heate d, shrunk) by a gro up tidal field after several gi- 
gayears (also see iMaver et al. II2OO1I), once the bound dark matter 
fraction drops below ~30% of its initial value. Even when only 
few percent of the stellar mass are lost, significant disk thickening 
(>50%) can occur for galaxies with low mass and/or high orbital 
ecce ntricities. 

iMcGee et all dioT l|) found that groups do not only contain 
a larger fraction of passive galaxies than the field, but that this 
fraction is larger for galaxies of lower mass and has grown over 
time, supporting the notion of continuous environmental influ- 
ence on dwarf galaxies in groups PI Furthermore, evidence for 
the releva nce of the l arge-scale environment is given by the find- 
ings of L ietzen et al.l (2012) that groups of equal richness have a 
higher fraction of elliptical galaxies when they reside in superclus- 
ter environments. Therefore, the fact that massive galaxy clusters 
contain, on avera ge, a larger fraction of d warf elliptical galaxies 
than groups (e.g. ITuUv & T rentham 2008) should not only be as- 
cribed to stronge r environmental effects inside the cluster (also see 
iBahe et alj|2013l) . Instead, the galaxy content of accreted groups 
had probably evolved further than that of groups in lower-density 
environments, which survived until today. 

When using our analysis of subhalo histories to interpret 
galaxy populations, we also need to consider that the baryonic 
structure at early epochs was different than in late-type galax- 
ies today. At 2: > 2, when part of the progenitors of today's 
dwarf ellipticals already experienced strong mass loss of their 
subhaloes (Fig. [8j, the di sk component may no t yet have had 
completed its formation dOovemato et al.l I2OIOI) . At redshifts 
2: >, 1, many star- for ming galaxies had a clumpy appearance 
dElmegreen et al.l |2004|) and likely represented an early phase 
of the spiral galaxy formation process (Boumaud et al. 2007|; 
lElmegreen et al.ll2009l) . These aspects would also be relevant for 
simulations of early interactions of gas-rich galaxi es that yield 
tidal dwarf galaxies dOabringhausen & KrouDall20T3l) . which may 



^ We note that this trend arises as a direct consequence of mapping mass 
loss to morphological type, since the anticorrelation of mass loss with clus- 
tercentric distance is stronger than with morphological type. 



We note t hat, for galaxies with luminosities equal or larger than 0.3L,, 
iBerrier et al.i ( l2009i) concluded that pre-processing is of minor importance. 
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constitute a re levant part of the observed dwarf galaxy population 
dKroupj 12012*). After their stellar structure had develop ed, dwarf 
galaxie s rarely experienced a major merger (De L ucia et aP 
l2012ah . Of the model galaxies in our SAM-V sample, only 3.1% 
have experienced at least one major merger since 2 = 1 (5.3% 
since z — 2)0 Major mergers thus only had a small contribu- 
tion to the structural appearance of dwarf galaxies in clusters today. 

Resulting from the interplay of various processes integrated 
over cosmic time, the present-day gala xy populations can be de- 
scribed by the classification scheme of Ivan den Berghl ( 1 19761) and 
iKormendv & Bended l l2012l) . This scheme groups lenticular and 
early-type dwarf galaxies into a parallel sequence to spiral and 
irregular galaxies. Given our analysis and the above considera- 
tions, today's early-type dwarf galaxies have not only experienced 
stronger subhalo mass loss than today's late-type galaxies of the 
same cluster, but have also resided in different environments for 
a large part of their lifetime. Therefore, the two parallel classi- 
fication sequences do not mean that the progenitor of a present- 
day early-type galaxy looked like a present-day late-type galaxy. 
Instead, the early-type sequence is most likely a consequence of 
stronger halo clustering and stronger influence of various envi- 
ronments over many g igayears — described as "history bias" by 
iDe Lucia et alj ( 1201 2bl) . The progenitors of the late-type sequence 
probably formed under different conditions, evolved at a different 
pace, and remained largely undisturbed by external influence. Late- 
type galaxies that are currently falling into massive clusters may 
thus not be good representations of the real progenitors of dwarf 
elliptical galaxies. 



6 CONCLUSIONS 

We investigated the history of the present-day dwarf galaxy popu- 
lation in clusters from the perspec tive of the Millenn ium-II ACDM 
cosmological simulation and the IGuo et al] | |201 lb semi-analytic 
galaxy model. In cases where a dark matter subhalo in the simu- 
lation did not lose a substantial amount of its mass over time, we 
can assume that the stars and gas were not affected noticeably by 
tidal forces, and consequently we identify these cases with late- 
type galaxies. Only in cases with strong subhalo mass loss can 
we assume that the baryonic configuration of the galaxy has also 
been affected significantly. These cases are identified with ellipti- 
cal galaxies, subsuming dwarf ellipticals and low-luminosity ellip- 
ticals. Assigning galaxy types to subhalo populations in this way 
leads to a morphology-distance relation very similar to what is ob- 
served in the Virgo and Coma clusters. 

We find that the median subhalo mass loss decreases more 
steeply with clustercentric distance than the observed elliptical 
fraction, and therefore also steeper than the increase of the late-type 
fraction. This can be resolved when assuming that ellipticals in the 
cluster core experienced even stronger environmental effects than 
ellipticals in the outskirts; the same applies to the intermediate-type 
and late-type galaxies. However, the properties of different model 
clusters and their projections show significant scatter, which could 
alleviate the differences. 

When subdividing them into early, intermediate and late types by the 
time spent in massive haloes, as in Fig.|5] the fractions ai'e 0.5% (2.1%) for 
early types, 2.4% (5.8%) for intemiediate types, and 5.9% (8.2%) for late 
types. These numbers increase only mildly when restricting our sample to 
the brightest dwarfs (-19.0 < Mr < -18.0 mag). 



The above statements also hold when using the time that sub- 
haloes have spent in massive clusters, instead of their mass loss. 
Similar concl usions have been drawn based on infal l time into a 
massive halo jPe Lucia et al.l2012bl : ISmith et al.l2012d ). While we 
argue that subhalo mass loss is more directly related to the tidal in- 
fluence on the baryonic galaxy, we also note that the cosmological 
simulation is limited in its particle resolution of subhaloes, and can 
therefore only provide approximate values. 

Our study shows that the majority of present-day dwarf ellipti- 
cals have already experienced strong mass loss of their subhaloes at 
high redshifts (z > 1). They have accordingly spent most of their 
lifetime in massive haloes (A/ lO^^MQ/h). This emphasizes 
the importance of environmental effects that acted in galaxy groups 
and "pre-processed" the galaxies before they entered the cluster. 
Dwarf ellipticals were thus not formed recently, but are likely a 
product of early and continuous environmental influence. We argue 
that this does not contradict small fractions of young stellar popu- 
lations in these galaxies. The central gravitational potential of the 
brighter early-type dwarfs is sufficiently deep to shield gas from 
being stripped. 

Over their lifetime, present-day late-type galaxies have expe- 
rienced an amount of environmental influence that is comparable to 
what the progenitors of dwarf elliptical galaxies had already expe- 
rienced at redshifts z > 2. In fact, there is no redshift at which the 
distributions of subhalo mass loss of today's late types and high- 
redshift progenitors of dwarf ellipticals agree. This reflects the fact 
that they evolved in different local and large-scale environments. 
Simulations aiming at reproducing the formation of dwarf elliptical 
galaxies would therefore need to take into account the environmen- 
tal characteristics of (proto-)clusters and groups at high redshift, as 
well as the fact that the progenitor galaxies themselves were at a 
much earlier stage of their evolution. 
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